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THE INFERENCE OF ATMOSPHERIC OZONE USING SATELLITE 


NADIR MEASUREMENTS IN THE 1042 cm" 1 BAND* 

By James M. Russell HI 
Langley Research Center 

and S. Roland Drayson 
University of Michigan 

SUMMARY 

Only a small amount of information is contained in radiances measured vertically 
from a satellite in the 1042 cm'l spectral region. Because of this low -information con- 
tent, it is not necessary to use an elaborate function to characterize the vertical ozone 
distribution. A function was used that described the ozone profile in terms of just three 
parameters — p^ 03 ^ m , p m , and H. These correspond, respectively, to the maximum 
ozone partial pressure, the atmospheric pressure at the altitude of the maximum partial 
pressure, and a width factor. An inversion procedure was designed to select values for 
these parameters so that the computed ozone spectrum at the top of the atmosphere 
matched the measured spectrum in a least-squares sense. It was shown that the calcula- 
tions should include the entire ozone band to obtain a maximum reduction of the effective 
random radiance noise. The atmosphere was divided into two broad layers: an upper 
layer containing pure ozone and a lower layer devoid of ozone containing the earth's sur- 
face and all the atmospheric water vapor and clouds. The problems connected with com- 
puting the water-vapor continuum and with considering clouds in the analysis were cir- 
cumvented by following this procedure. The top of the lower layer was taken as the 
boundary under the satellite, and the effective temperature of the layer was used to com- 
pute the lower boundary emission in the equation of radiative transfer. The effective 
temperature could be found to be within 0.5 K or less by considering water -vapor line 
absorption in the calculations. 

The least-squares solution for the three ozone parameters was expanded in terms 
of the eigenvectors and eigenvalues of a solution matrix. A subsequent analysis using 
synthetic clear -sky radiances showed that, under favorable conditions, measurements in 
the ozone band contain two pieces of independent information: the total ozone u, and the 
altitude h m (or pressure p m ), at the level of the maximum ozone partial pressure. An 

The basic information presented herein was included in a thesis by the first author 
and was offered in partial fulfillment of the requirements for the degree of Doctor of 
Philosophy, University of Michigan, Ann Arbor, Michigan, July 1970. 



error study revealed that errors in u are affected most by random radiance errors, 
lower boundary temperature errors, and ozone absorption-line intensity errors. Errors 
in h m are affected most by the former two errors and by temperature -profile bias 
errors. When all error sources were considered simultaneously, the results showed that 
it should ultimately be possible to determine the total ozone to within 10 percent or less 
and to determine the altitude of the maximum partial pressure in the ozone profile to 
within 1.5 km when the root -mean -square radiance noise level is 1 percent or less. 

These latter calculations were repeated for varying degrees of cloudiness in the tropo- 
sphere. The data showed that the presence of clouds did not seriously affect results as 
long as there was some contrast between the ozone spectrum and the effective lower 
boundary emission spectrum. 

Next, the inversion technique was applied to balloon data measured over Palestine, 
Texas, and to the Nimbus HI satellite data measured over Grand Turk Air Force Base in 
the Bahamas. The Palestine data contained very low radiance noise, and the correspond- 
ing inversion gave a value for u which was within about 3 percent of the value estimated 
from a nearby Dobson measurement. The value for h m could not be accurately esti- 
mated, but the result obtained from the inversion was within the possible range of values 
indicated by past measurements made at the same general location and time of year. The 
Nimbus HI data were very noisy by comparison to the balloon measurements, and this 
noise was reflected in the inversion. The inferred values for u and h m were off by 
76 percent and 6.3 km, respectively, when compared to values determined by a simulta- 
neous ozonesonde observation. These errors are consistent with the predicted values of 
96 percent and 7.7 km, respectively, as determined from the error-study results and the 
estimated random radiance noise in the Nimbus III data. 

The application of the inversion technique to the Palestine and Nimbus data proved 
the validity of the inversion method and demonstrated the usefulness of the error study in 
predicting errors. The calculations also illustrated a critical need for further work to 
improve the accuracies of the absorption-line intensities. The basic difference between 
the Palestine and Nimbus data was the radiance noise level; the noise for the Palestine 
radiances was low, whereas the noise for the Nimbus data was comparatively high. It is 
believed that this is the underlying reason for the drastic differences in the inversion 
results for the two sets of data. If the absorption-line intensities had been known with 
greater precision, a much wider spectral interval could have been used in the calcula- 
tions, thereby causing a significant reduction of the effective noise level in the Nimbus HI 
data and a substantial improvement in the results. These results indicate that, unless 
absorption errors are decreased by empirical adjustments or other methods, calculations 
will have to be limited to a narrow spectral interval and, consequently, to low-noise data 
(less than 1 percent). 
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INTRODUCTION 


Ozone is only a minor atmospheric constituent; yet, it has held the attention of sci- 
entists ever since it was first suggested by the chemist Schonbein in 1840. The first con- 
clusive evidence of an ozone layer in the atmosphere was deduced from ultraviolet absorp- 
tion measurements made by Fabry and Buisson (ref. 1). Later, their conclusions were 
given theoretical support by the British geophysicist Chapman (ref. 2) who first explained 
the presence of ozone in the upper atmosphere by using photochemical principles. At 
about the same time, Dobson (ref. 3) published his paper describing the development of 
Fabry and Buisson's technique into a convenient method for routinely measuring the total 
ozone above an observer. 1 

After the work by Dobson and his coworkers the real meteorological significance of 
atmospheric ozone began to emerge. The systematic measurements of total ozone made 
possible by the Dobson instrument revealed a strong correlation between ozone amount 
and day-to-day weather variations. Later, with the advent of rocketsondes and high- 
altitude balloons, instruments were developed to measure, in situ, the vertical distribu- 
tion of ozone as well as temperature and wind-field profiles. These measurements 
showed definite correlations between significant changes in the ozone partial pressure and 
changes in the temperature lapse rate and the wind-field frequency. (See, for example, 
ref. 4, or, for an extensive bibliography, ref. 5.) Other measurements showed an unques- 
tionable dependence of total-ozone content on the position of the subtropical jet stream 
relative to the location of the ozone observations (refs. 6 and 7). Measurements by 
Godson (ref. 8) revealed correlations on a finer scale. He compared the 10-day running 
mean of 100-mb temperature with the total ozone for data taken over a 9-month period at 
Edmonton, Canada. His results showed that the average daily temperature variation was 
almost exactly correlated with the average daily variation in ozone amount. 

As measurements of the vertical distribution of ozone accumulated and the photo- 
chemical theory was developed, important discrepancies were noticed between the distri- 
butions observed and the profiles predicted based on photochemistry. These discrepan- 
cies gave the first indication of the significant interrelationship between ozone content and 
atmospheric motions. Wulf (ref. 9) was apparently the first to propose a general circula- 
tion model which explained the ozone spring polar maximum. After his work, other 
papers were published concerning the relation between ozone and atmospheric motions, 
and the role of ozone in this regard became widely accepted. (See, for example, refs. 10 
to 13.) 

*This is the total (integrated) ozone amount in an atmospheric column of unit cross 
section which extends from the ground to the top of the atmosphere. The phrases "total 
ozone" and "ozone amount" will be used synonymously throughout this report. 
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The complex interplay between ozone and other atmospheric properties, especially 
atmospheric motions, caused widespread interest in global measurements of both total 
ozone and the vertical ozone distribution. In 1965, the National Academy of Sciences pub- 
lished a report outlining goals for an international observation program (ref. 14). The 
report cited several areas of research which would be advanced by global ozone measure- 
ments. These areas include the study of atmospheric dynamics in the region of the upper 
troposphere and lower stratosphere, the analysis and study of the circulation patterns of 
the Northern and Southern Hemispheres and the exchange processes that take place 
between them, the study of the mechanics of the tropospheric-stratospheric exchange 
process, the study of stratospheric photochemistry, and the study of atmospheric phenom- 
ena such as the "explosive warming" of the winter-spring stratosphere and the 26-month 
oscillation in the stratospheric wind and temperature. 

The general interest in global- scale measurements of atmospheric ozone resulted 
in the conception of a number of satellite measurement techniques. These methods 
include the measurement of backscattered solar ultraviolet (UV) light (refs. 15 to 20), the 
measurement of visible and near UV solar radiation passing tangentially through the 
atmosphere (refs. 21 to 23), the measurement of the UV spectrum of a star during occul- 
tation from a satellite (ref. 24), and the use of laser transmission from the ground to a 
satellite or vice versa (ref. 5). The use of satellite nadir measurements in an infrared 
band (1042 cm -1 ) of ozone has also been suggested (ref. 25). This method is attractive 
since it does not contain many of the problems of the other techniques. In the infrared 
technique, the measurements are not restricted to specific geographic areas, and the 
method is applicable on both the day and night side of the planet. Also, since scattering 
is unimportant in the infrared spectral region, the analysis is comparatively simple and 
the multiple-scattering problem of the UV backscatter method does not exist. 

Studies dealing with the determination of atmospheric ozone amounts from satellite 
nadir measurements in the infrared have been reported by Sekihara and Walshaw (ref. 26) 
and Prabhakara (ref. 27). A more recent paper by Prabhakara et al. (ref. 28) discusses 
inversion results for radiance measurements made with an Infrared Interferometer Spec- 
trometer (IRIS) from on-board the Nimbus III satellite. In the latter report, the atmo- 
spheric ozone profile is represented by using empirical orthogonal functions. The atmo- 
spheric transmission due to ozone is computed by using a statistical band model, and the 
absorption-line parameters used in the computations are determined by trial and error 
with the help of an ozone sounding and Nimbus III radiance measurements. The water- 
vapor continuum is considered, but absorption by water-vapor line centers is not. 

The present report presents the results of a detailed study of the nadir experiment. 
A different approach from that of Prabhakara et al. (ref. 28) is developed which reveals 
explicitly how many pieces of information are available from the measurements. The 
statistical band model is used to compute transmittances for an error study; but when the 
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inversion procedure is applied to balloon and satellite data, the line-by-line method is 
employed. No empirical corrections are imposed on the absorption-line parameters 
used in these calculations. A method is also presented to account for the presence of 
water-vapor line absorption in the ozone band. 
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SYMBOLS 

Planck radiance, watts/meter-steradian 
velocity of light in a vacuum, centimeters/second 

mean spacing of absorption lines in an interval Av wide, per centimeter 
ground state energy per centimeter 

ratio of perturbed to unperturbed line strength, dimensionless 

absorption-line profile function, centimeters 

acceleration of gravity, centimeters/second 2 

width factor in Green's ozone profile function, dimensionless 

altitude of maximum ozone partial pressure, kilometers 

total number of atmospheric layers, dimensionless 

Boltzmann constant, ergs/molecule-kelvin 

Planck radiance for effective lower boundary temperature, 
watts / mete r - ster adian 

spectral radiance, watts/meter-steradian 

effective length of a homogeneous absorption path, centimeters 
mass per molecule, grams 

total number of absorption lines in a band, dimensionless 
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total number of absorption lines in an interval Av wide, dimensionless 
number density at standard temperature and pressure, per cubic centimeter 
probability distribution function for line intensity, dimensionless 
probability distribution function for line positions, dimensionless 
atmospheric pressure, millibars 

atmospheric pressure at level of maximum ozone partial pressure, millibars 
ozone partial pressure, micromillibars 

maximum ozone partial pressure in Green's ozone profile function, 
micromillibars 

mass mixing ratio of an absorbing gas, dimensionless 
intensity of ith absorption line per atmosphere-centimeter 2 
atmospheric temperature, kelvins 
variable of integration, dimensionless 

integrated ozone mass in an atmospheric column of unit cross section, 
atmosphere-centimeter 

ith eigenvector, unit varies 

a constant, dimensionless 

Doppler half-width, per centimeter 

Lorentz half-width, per centimeter 

mean absorption-line intensity for interval Av per atmosphere-centimeter 2 

unperturbed line strength in ground-state representation per 
atmosphere-centimeter 2 
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6g(i>) emissivity of lower boundary under a satellite, dimensionless 

k(P,p) absorption coefficient, per centimeter 

X wavelength, centimeters 

Ai ith eigenvalue, unit varies 

v wave number, 1/X, centimeter -1 

a standard deviation of a random variable, unit varies 

r(D,p) transmittance, dimensionless 

0( v) radiometer instrument function, dimensionless 

ANALYSIS 

Introduction 

A spectrum of the upwelling thermal energy measured from a balloon above 
Palestine, Texas, on May 8, 1966, is shown in figure 1 (refs. 25 and 29). These data 
were obtained with a breadboard version of an IRIS in a test conducted by the High Alti- 
tude Engineering Laboratory of the University of Michigan for the Goddard Space Flight 



Figure 1.- Atmospheric emission spectrum acquired with a breadboard version of an ERIS 
instrument in a balloon flight at an altitude of 31 km on May 8, 1966, from 
Palestine, Texas (after Conrath, ref. 25). Smooth curves are plots of blackbody 
function at the indicated temperatures . 
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Center. The spectrum shows marked structure due to absorption by carbon dioxide, 
ozone, and water vapor. The magnitude of the energy at any given wave number in such 
a spectrum is a function of several variables. It depends on the temperature of the lower 
boundary (if the atmosphere is not opaque in the wave-number region of interest), on the 
temperature profile in the atmosphere, and on the amount and vertical distribution of 
gases that emit and absorb energy in the spectral pass band of the instrument. The 
upwelling energy also depends on the infrared reflectance properties of the lower bound- 
ary, but this dependence can usually be neglected. The fundamental relationship between 
the outgoing energy and properties of the earth -atmosphere system is the equation of 
radiative transfer (ref. 30). The form of the equation used in this paper is based on the 
following assumptions: 

(a) The atmosphere is plane -parallel. 

(b) Scattering at the wavelengths under consideration can be neglected. 

(c) Atmospheric refraction can be neglected. 

(d) The absorbing gases are in local thermodynamic equilibrium up to an altitude of 
60 or 70 km (ref. 31). 

If it is further assumed that the satellite instrument receives energy in a narrow cone in 
the vertical so that everywhere in the cone sec 6 ~ 1 (where 0 is the zenith angle), 
the equation of radiative transfer can be written in integral form as 



dT(5> 0 ,p) 

d(log p) 


d(log p) (1) 


where L^ 0 0) is the radiance at the top of the atmosphere for a given wave number v Q , 
T(p) is the atmospheric temperature, B^ 0 ,T(p)j is the Planck radiance, eg(^ 0 ) ^ 

Pg are the emissivity and pressure, respectively, at the lower boundary, and r(v 0 ,p) is 
the transmittance of the atmosphere from the satellite down to an arbitrary height where 
the atmospheric pressure is p. The bars above the symbols denote that these quantities 
are average values and the o subscript denotes the wave number at the center of the 
averaging interval. The averaging is performed according to some function which ade- 
quately represents the frequency response of the instrument. 


The first term on the right in equation (1) represents the contribution to the radi- 
ance at the top of the atmosphere due to emission by the lower boundary (e.g., the earth's 
surface). The second term accounts for the contribution made by atmospheric emissions. 
Any descending energy which is reflected at the lower boundary has been neglected. 

Atmospheric ozone information enters equation (1) through the transmittance 
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Obviously, ozone data cannot be retrieved by using equation (1) unless the atmospheric 
temperature profile is known. It is assumed in this study that temperature data will be 
available from any one of the numerous satellite remote -sensing techniques which have 
been proposed (e.g., refs. 32 to 34). Before considering the problem of inverting equa- 
tion (1) to obtain ozone information, there are two important topics which must be dis- 
cussed: (1) the method used to compute the atmospheric transmittance, and (2) the 
method used to compute the lower boundary emission. 


Atmospheric Transmittances 

The transmission of the atmosphere in the vicinity of 1042 cm - * is influenced by 
water vapor and carbon dioxide as well as ozone. The water-vapor continuum absorption 
must be considered as well as absorption by water-vapor lines superimposed on the con- 
tinuum. However, carbon-dioxide absorption is weak and it will be neglected in the cur- 
rent calculations. The estimated upper limit to the transmittance error caused by this 
assumption is about 5 percent based on calculations using the weak-line approximation. 

The actual error is probably only about 2 to 3 percent. Carbon-dioxide effects should not 
be neglected in future computations once the accuracies of the ozone transmittances are 
improved. 

Water-vapor line absorption can be accurately computed if the vertical water-vapor 
distribution is known; and even if the vertical profile is not known, an adequate treatment 
is possible. However, the water-vapor continuum spectrum is difficult to model analyti- 
cally since it arises from absorption in the wings of many lines located in distant-rotation 
and vibration-rotation bands. There are presently no totally acceptable theories avail- 
able concerning this problem. Consequently, the influence of water vapor on the atmo- 
spheric transmittance in the ozone band will not be computed directly in this study. 

Instead, the effect will be treated as part of the problem of determining the lower bound- 
ary emission. The way this is done is discussed in the next subsection. Thus, the prob- 
lem of computing the transmittance of the atmosphere is, in a sense, reduced to one of 
computing the transmittance due to ozone alone. Two methods of accomplishing this are 
described in the subsequent paragraphs. 

The atmospheric transmittance at a wave number v from the top of the atmosphere 
down to a height where the atmospheric pressure is p is given by 


T(D,p) = exp - 



(2) 


where q(p) is the mass-mixing ratio of the absorbing gas, g(p) is the acceleration of 
gravity, and /c^p) is the absorption coefficient. Assume now that the atmosphere is 
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divided into small incremental slabs in the vertical and consider the transmission through 
only one of these layers. If each layer is considered to be homogeneous, the transmit- 
tance for the jth layer is 



(3) 


where pj is the mean pressure for the layer and Uj is the integrated mass across the 
layer given by 


_ PP2 q(p) 
j 1. g(p) P 


Pi 


(4) 


The pressures at the boundaries of the layer are given by p^ and P 2 - Equation (3) is 
not an exact expression since, in reality, an atmospheric slab is inhomogeneous and the 
absorption coefficient varies with pressure across the layer. However, errors which 
arise in using equation (3) can be made as small as desired by using very thin slabs. 
Fundamentally, the problem of computing the transmittance of the atmosphere is one of 
calculating the absorption coefficient for each of the j layers in the atmosphere. This 
is not a simple problem since the absorption coefficient at any given wave number v is 
the sum of the coefficients at v due to each of a large number of absorption lines in a 
band, and each of these lines varies with pressure and temperature. 

Consider first the absorption coefficient at v due to a single line centered at a 
wave number ^ in the band. For this ith line and the jth atmospheric layer, 

«i(^,Pj) = % f(^i) (5) 

where Si is the line intensity or strength and is the profile function which 

describes the line shape. The profile function is normalized so that 

jo” - 1 < 6) 

and 

^ = j'o" “i < 7 > 

The absorption coefficient varies with altitude because of changes in both the line 
strength and the profile function. The line strength contributes to the variation through a 
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temperature dependence. The relationship used by Clough and Kneizys (ref. 35) for ozone 
is given by 


Si(Tj) = XnFi^j 


5/2 


exp 


fo 

kT j 


■J 


( 8 ) 


where X is a constant, y j is the unperturbed line strength in the ground-state repre- 
sentation, Fj is the ratio of perturbed to unperturbed line strength, E 0 is the ground- 
state energy, and k is the Boltzmann constant. The profile function varies with altitude 
because of both a temperature and a pressure dependence. The variation in the profile 
function for the range of physical conditions considered here can be described by two 
limiting cases: the Lorentz and Doppler line shapes. A third case, the natural-line 
shape, can be neglected since the width of the natural broadened line is much smaller than 
the width for either of the other two cases. 

The Lorentz line shape is caused by molecular collisions and is described by the 
function 



where is the Lorentz half-width at half-maximum. The half-width depends on 

pressure and temperature according to 



( 10 ) 


where agi is the half-width at temperature T 0 and pressure p Q . 

The Doppler line shape arises because of thermal motions of the molecules. In 
this case, 


where 





x = 




( 11 ) 
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and where a 


is the Doppler half-width at half-maximum which is expressed by 


Di 
<*Di 

The quantity m is the mass per molecule in grams. 

The transition from the region of the atmosphere where one limiting shape is 
important to the region where the other dominates is not abrupt. Therefore, a range of 
pressures is encountered where both thermal and collisional effects are important in 
determining the absorption-line profile. Under these conditions, the correct line shape 
is given by a convolution of the Doppler and Lorentz shapes. This shape, called the Voigt 
function, is represented by 


V 2kT j 

— In 2 J 

c \ m . 


, 1/2 


( 12 ) 



So far, the main concern has been with a single absorbing line and it has been shown 
how the absorption coefficient for this line varies as a function of the atmospheric param- 
eters, pressure, and temperature. To compute the integrated effect of all lines in a band 
on the absorption at a given wave number P, it is necessary to sum the values of absorp- 
tion coefficients at v due to each absorption line. Thus, for N lines, the monochro- 
matic absorption coefficient is 


N 

K'^i) = 2 ) 

i=l 


(14) 


The monochromatic transmittance is computed by inserting equation (14) into equation (2); 
and then the average value is calculated according to the following equation: 


t (5/p) 



T(v,p)cp{l>)dv 



(15) 
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where and Pg are the bounds of the spectral interval spanned by the instrument 
function Equation (14) is a lengthy calculation because of the large number of 

lines normally found in an absorption band. The equation is also complicated because of 
the requirement that the computations be done for all j layers in the atmosphere, in 
which case the line strength and shape variability must be included. These difficulties 
prompted the development of the band-model concept to compute transmittances. In 
this approach, the actual distributions of absorption-line positions and intensities are 
simulated; and, instead of computing equation (14) and the average transmittance from 
equation (15) directly, the calculations are approximated by an analytical integration 
using the assumed distributions of the line parameters. Even though band models are 
still frequently used today, it is now possible, because of the existence of high-speed 
computers, to perform the calculations directly by the direct-integration (or line-by-line) 
method (refs. 36 and 37). The advantages and disadvantages of each of these approaches 
were discussed by Drayson (ref. 38) and will not be repeated here. Both procedures 
were used in this present work and a brief description of each method follows. 

The direct-integration method is numerically the most accurate of the two tech- 
niques. The method predicts the absorption almost exactly if the absorption-line param- 
eters are known precisely. In the program developed by Drayson and Young (ref. 37) the 
monochromatic transmittance is computed at various wave numbers in the absorption 
band. Since the values are computed monochromatically, the transmittance from the top 
of the atmosphere down to a height where the pressure is p can be computed by multi- 
plying the transmittances for each of K atmospheric layers. Accordingly, 


t(«>,P) 



(16) 


where the symbol n denotes multiplication. These monochromatic transmittances are 
then averaged over 0.1 cm _ l intervals at the bottom of each Kth layer. With this tech- 
nique, the atmosphere can be divided into a series of thin homogeneous layers, each of 
which has values of pressure and temperature equal to the average values for the layer. 
Since each transmittance value is computed for a thin atmospheric layer, it is relatively 
easy to include the appropriate absorption-line shape and intensity for that region of the 
atmosphere in the calculations. The Drayson and Young program includes all three 
absorption-line broadening regimes, it allows for the variation of line half-width with wave 
number, and it accounts for the variation of line intensity with temperature. The chief 
drawback to using the method is that it is very time consuming, even when modern, high- 
speed computers are used. Consequently, the technique was used here only when a high 
degree of computational accuracy was desired, such as in calculating the effect of line- 
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strength temperature dependence on transmittance and in applying the inversion method 
developed in this paper to the Nimbus HI and Palestine, Texas, flight data. The error- 
study transmittances were calculated by using the random exponential band model 
(refs. 39 and 40). A description of this model follows. 

The major gain in using the band-model technique rather than the line-by-line 
method is a great reduction in the time required for a calculation. The line-parameter 
distribution functions used with the most success in the ozone band are a random distri- 
bution of line positions and an exponential distribution of line intensities. (See, for exam- 
ple, ref. 41.) The probability distribution function for the line positions is P(P) = -^r 

Av 

where Av is the width of the spectral interval, and the function for line intensities is 

P(S) =i-exp^|-^ (17) 

where is the mean absorption-line intensity for the interval Av. With these distri- 
butions, the average transmittance at a wave number v through a single atmospheric 
layer is (ref. 42) given as 


^,Pj) = exp(-w) (18) 


The argument of the exponential is computed as follows: 



(19) 


where d = — is the mean line spacing and N' is the number of lines in the interval 
N’ 

Av. Since the transmittance computed with the band model is an average value for the 
spectral interval, the transmittance through K layers cannot be computed simply by 
multiplying the values for each jth layer (ref. 42, p. 123). Instead, the transmittance 
must be calculated for an equivalent homogeneous slab defined by the K layers. The 
reduction of an arbitrary path in the atmosphere to an equivalent homogeneous path is 
usually done by means of the Curtis-Godson approximation (e.g., ref. 43). This approxi- 
mation is not very accurate for atmospheric ozone (refs. 44 and 45), but the accuracy is 
good enough that the method can still be used in the error-study calculations discussed 
later. There are two other approximations which will be used in the band-model calcula- 
tions. It is assumed that purely collisional line broadening is valid throughout the range 
of altitudes considered and that the line intensities do not vary with temperature. The 
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assumption of collisional broadening causes a small transmittance error above an alti- 
tude of approximately 30 km since in this region the Voigt profile function is important 
(ref. 31). However, the estimated maximum error is only about 2 percent according to 
the correction table published by Gille and Ellingson (ref. 46). The effect of neglecting 
the line -intensity temperature dependence is discussed in a subsequent section. These 
approximations cannot be used in precision calculations but they should have little effect 
upon the error-study results. Goody (ref. 42) integrated equation (19) for a Lorentz pro- 
file function to give 



= exp 



^ k 
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~\ 1/2 

Jl + 2 

" H 1 

?k ] 

1 

2 ™Lo (I> ) 

i J 


( 20 ) 


where and u^ are the effective length and absorber amount, respectively, of the 
homogeneous path and (P) is the Lorentz half-width at the wave number u and at 
standard temperature and pressure. The subscript k refers to the effective value for 
a homogeneous slab defined by K layers. The methods for computing the various quan- 
tities in equation (20) are discussed in appendix A. A spectral interval of 5 cm -1 was 
used in all calculations with equation (20) to simulate the resolution of an IRIS. 

All transmittance calculations were based on ozone absorption-line strength and 
position data published by Clough and Kneizys (refs. 35 and 47). A good indication of 
errors that may exist in these data can be obtained by comparing theoretical and experi- 
mental absorption spectra for a homogeneous path. S. R. Drayson and C. Young of the 
University of Michigan have done this by using their line-by-line technique and experi- 
mental measurements made by McCaa and Shaw (ref. 48). The experimental values are 
believed to be accurate to within 10 percent or less. A comparison of the computed and 
experimental spectra for a set of conditions similar to those encountered in the atmo- 
sphere is shown in figure 2. The theory predicts too little absorption on the low wave- 
number side of the band and too much absorption on the high wave-number side. Also, 
the discrepancies are much greater than the experimental inaccuracies. It appears, 
therefore, that there are large errors in the absorption-line parameters. Nevertheless, 
no attempt is made in this study to correct the values of these parameters either empir- 
ically or by other means. Empirical adjustments can be a useful means of correcting the 
absorption parameters as shown by Prabhakara et al. (ref. 28). However, corrections 
made for atmospheric conditions prevailing above one region of the globe may not neces-' 
sarily be appropriate corrections for a different region. Ultimately, therefore, it 
appears that the most satisfactory approach is to use absorption parameters based only 
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Figure 2.- Comparison of experimental and calculated absorption 
in the 1042 cm -1 ozone band. 


on a knowledge of the physics of the ozone molecule as is done in this report. The effects 
of absorption-parameter errors on results are discussed in the error-study section of 
this paper. 


Lower Boundary Emission 

Background emission from the underlying boundary can significantly alter the ozone 
spectrum measured at the top of the atmosphere. Sometimes, in high-latitude areas of 
the globe, the combined effect due to atmospheric and lower boundary emissions can 
create the condition where there is little or no contrast between the background emission 
spectrum and the ozone spectrum. Under these conditions, it may be difficult or impos- 
sible to determine ozone information from satellite measurements. These effects are 
discussed in appendix B. The following subsection describes the technique used in this 
study to compute the lower boundary emission. 

The lower boundary emission is clearly represented by the first term on the right- 
hand side of equation (1). However, the evaluation of this term is not straightforward. 

For example, if the underlying boundary is the earth's surface, it is difficult to determine 
the surface conditions which should be used. The field of view of a satellite instrument 
covers a large area; and since both the surface temperature and surface emissivity vary 
with terrain, appropriate average values must be found for these quantities. If the sky 
contains clouds, the emissivity of the clouds must be considered and the emissions from 
the clouds and surface must be weighted in some way according to the percentage of cloud 
cover. The problem is further compounded by the presence of the water-vapor continuum 
absorption which is difficult to compute. This absorption attenuates the energy emitted 
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by the earth's surface before it reaches the top of the atmosphere. But these difficulties 
can be circumvented to a large extent by using the fact that the ozone partial pressure is 
small in the troposphere. This property of atmospheric ozone makes it possible to 
divide the atmosphere into two broad layers: an upper layer containing pure ozone and a 
lower layer devoid of ozone containing the earth's surface and all the atmospheric water 
vapor and clouds. This approach greatly simplifies the problem as will be seen in the 
subsequent paragraphs. 

If the spectral interval is sufficiently wide so that there is no correlation between 
O3 and H2O absorption, then the transmission through a water-vapor ozone atmosphere 
is given by (ref. 42, p. 123) 


t(v, P ) = t(h 2 0,P,p) t(o 3 ,v, P ) (21) 

If, for clarity, the arguments of the variables in equation (1) are dropped and the sub- 
script g is used to denote the lower boundary and 0 to denote the top of the atmo- 
sphere, the equation can be written 


L 0 = egBgr g (H 2 0) T g (o 3 ) + ^ B df^HgO) r(o 3 ) 

Trr L_ _ 


( 22 ) 


where the relation of equation (21) is used. Expanding equation (22) gives 
L 0 = £g B g r g (H 2 0) r g (o 3 ) + B?(03)df(H 2 0) ♦ £ , BT^drfos) 


(23) 


Consider now figure 3 which shows typical water vapor, ozone, and temperature distribu- 
tions and which illustrates the two-layer concept. Since the water-vapor content is very 
small in the stratosphere and the 0 3 content is small in the troposphere, t^ 0 3 ^ - Constant 

in the lower layer and t^H 2 6) - 1 (in the spectral vicinity of 1042 cm - *) in the upper 

region. For these conditions ^(03) ~ ^(03) where ^ r (o^ is the transmittance from 
the top of the atmosphere down to the bottom of the ozone layer. Therefore, equation (22) 
can be written 


Ln — 


e g B g T| 


■sM + £‘ 
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;(h 2 o) 


B d 


5(H 2 o) 


( 03 ) + l! r ( 03 ) B<if ( O 3) < 24 > 
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Figure 3.- The two-layer concept. 


where the bracketed portion of the first term on the right-hand side is the emission from 
the lower layer. If the water -vapor line -absorption spectrum is excluded for the present 
and the emissivity (e^j is considered to be constant in the spectral regions immediately 
adjacent to and including the ozone band, the spectrum at the top of the atmosphere can 
be viewed as consisting of the Planck radiance of the lower layer at some effective tem- 
perature T e ff and the superimposed absorption spectrum of a pure ozone layer. The 
temperature T e ff will depend on the surface, water vapor, and cloud conditions. The 
value of T e ff is determined by selecting a value which gives a match between the Planck 
radiance and the measured spectrum on either side of the ozone band. 


Since the measurements contain noise, T e ff is found by using the method of least 
squares fitting a Planck function through the maximum radiance points of the measured 
spectrum. This is illustrated by the dashed line in the bottom spectrum of figure 4. The 
spectral structure in figure 4 is due to ozone absorption and line absorption by water 
vapor. It is not obvious that water -vapor line structure exists in the part of the spectrum 
where ozone absorbs. This is because the effects of the lines are reduced in the center 
of the band by absorption of the upwelling lower boundary energy in the ozone layer, and 
they are hidden in other spectral regions because of the slope of the ozone spectrum. 
Nevertheless, line structure is present throughout the band as part of the background 
emission. This is illustrated in the lower spectrum of figure 4 by the dotted water-vapor 
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LINE LINE 



Figure 4.- Sketch of typical ozone and water -vapor absorption spectra. 

.lines... ^The^magnitudes of these lines are not drawn to scale relative to the ozone spec- 
trum, and the positions do not necessarily coincide with the positions in an actual water- 
vapor spectrum.) The effect of water-vapor lines is considered by computing a line- 
absorption spectrum for tropospheric water vapor by using the direct-integration 
technique. The top spectrum in figure 4 illustrates this computation. Preferably, the 
calculations would be done by using an actual sounding of the water-vapor altitude distri- 
bution. However, if there is no sounding available, the mean climatological profile can 
be used. After computing the water-vapor line spectrum, the water-vapor absorption at 
any point in the measured spectrum can be predicted based on the relative absorption in 
the computed water-vapor line spectrum. For example, suppose it is desired to compute 
the background emission at the position of the line A in the measured spectrum. 

The ratio of the absorption by line A to that by line B in the computed line spectrum 
is a/b. Therefore, the lower boundary emission B^ at the position of water-vapor 
line A is 
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(25) 


BA = l a - ^ (l b - B B ) 


where La and L B are the Planck radiances at the positions of lines A and B for 
T = T e ff, and B B is the measured upwelling radiance at the position of line B. This 
method was used with the Nimbus and Palestine, Texas, flight data. The direct-integration 
method and the H2O line parameters of Benedict and Calfee (ref. 49) were used to compute 
the line absorption. Some calculations were done for water-vapor absorption centers 
located outside the spectral region of ozone absorption. It was possible with this proce- 
dure to evaluate the accuracy of lower boundary emission values determined by this 
method. The calculations showed that the background emission could be computed to 
within 5 x 10“^ watt/m-sr or less. This corresponds to an effective background- 
temperature error of =0.5 K in the wings of the ozone band and a smaller error toward 
the center. The background emission computed in this manner at various wave numbers 


can be used to replace quantity 



in the first term on the right-hand 


side of equation (1). When this is done, the transmittances in that equation correspond to 
the transmittances of the atmosphere due to ozone alone. 


The Inversion Method 

The problem to be solved is the following: given the temperature variation with 
altitude, find an atmospheric ozone profile which gives a match between the radiance 
l c(^oj°) which is calculated with equation (1) and the radiance L m ^ o ,0) which is mea- 
sured. Equation (1) behaves like a Fredholm integral equation of the first kind (ref. 32) 
Therefore, an infinite number of ozone profiles can be found which will provide a radiance 
match to within any given error limits. This family of results will include a multitude of 
solutions which bear no physical resemblance to a realistic ozone profile. Therefore, 
one is faced with the problem of filtering out those solutions which are mathematically 
valid but which are physically unrealistic. This is normally accomplished by preselect- 
ing a form of the solution which is known to be physically valid based on other informa- 
tion, that is, climatology. (For example, in the temperature-inversion problem discussed 
by Wark and Fleming (ref. 32), the temperature solution is expressed in terms of empiri- 
cal orthogonal functions which are determined from a large number of temperature sound- 
ings in a given area of the globe.) (See ref. 50.) 

A function is needed which adequately represents the general shape of the vertical 
ozone distribution described in any one of a number of references (e.g., refs. 51 to 54). 

It would be desirable to include enough parameters in the function to define the dominant 
features, that is, the magnitudes and altitudes (or pressure levels) of the three major 
peaks in the profile. However, because of the low -information content of the measure- 
ments (ref. 55), the most that can be expected from an inversion is data that define the 


20 



gross features of atmospheric ozone, that is, properties such as the total ozone, the alti- 
tude where the main peak in the profile occurs, and the value of the maximum partial 
pressure at the main peak. A function which characterizes the profile by just a few 
parameters is all that is required, since greater detail is not available from the mea- 
surements. Green (ref. 56) used a relationship that is well suited for this purpose. It is 
expressed in terms of the value of the maximum partial pressure p( 03 ) m > the atmo- 
spheric pressure p m at the altitude where the maximum partial pressure occurs, and 
a width factor H. The function is given by the following equation: 


where 


p(o 3 ) = 4 p(o 3 ) m 


(26) 


(i + e*y 


x = 


l°g(p/Pm) 


H 


and where p(o 3 ) is the ozone partial pressure. The function is used in this study to 
describe the ozone profile in a range of altitudes starting near the tropopause and extend- 
ing upward. The partial pressure in the troposphere is considered to be constant. A 
plot of this function is shown in a subsequent figure. 

The goal of the inversion method is to select values for the three parameters 
P(o 3 )m’ Pm’ 311(1 H which give a match between L c ^ o ,0^ and L m (P 0 ,(j) according 
to a criterion which is discussed subsequently. The initial step in the solution procedure 
is to make a linear approximation for the relationship between L m (P 0 ,d), L c (v o ,0), and 
the three parameters as follows: 


— /_ \ — /_ \ 9L (v,0) 
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(27) 


m 


This is a Taylor's series expansion with the higher order terms truncated. To determine 
three parameters, at least three equations must be used. The required three equations 
can be written from equation (27) to correspond to measured radiances at three wave 
numbers. Dropping the parentheses and using a numerical subscript to denote the first, 
second, and third wave-number centers gives 
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(28) 


An iterative procedure is used to determine the final values of the three parameters 
p(c>3)m5 P m > an< ^ H. After making initial guesses for the parameters, the quantities 

L c ^ are computed by using equations (26), (3) or (20), and (1). The partial derivatives 
are computed in a similar manner by letting one parameter at a time change by a small 
amount in order to obtain the differentials. Upon solving for Ap^O^^, Ap m , and AH, 
these quantities are added to the initial guesses and the process is repeated until some 
convergence criteria are met. 

The system of equations (28) could be solved directly without further manipulations. 

However, since the measured radiances will contain noise, it is desirable to use a least- 

squares fitting procedure. It is also desirable to use an overdetermined solution method; 

that is, use many more equations than are needed to determine the three parameters. By 

using a large number of measurements (preferably measurements covering the entire 

ozone band) and by using a least-squares solution, the standard deviation of the radiance 

o 

noise can, in effect, be reduced. So, if M measurements are used, instead of having 
just three equations, the system of equations (28) will have M equations. This system 
of M equations can be expressed in matrix form as follows: 

AL = Q R (29) 

where AL is an M x 1 matrix defined by the differences (^L m j, - L C i), Q is the 
M x 3 matrix of partial derivitives, and R is the 3x1 matrix defined by the pertur- 
bation terms. The least-squares solution to equation (29) is easily shown to be 

3 Currently, it is not possible to use the entire band because of the magnitudes and 
differences in sign of absorption errors in various parts of the band. (Refer to fig. 2.) 
This emphasizes the importance of further work to improve the accuracies of the 
absorption-line parameters. 
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(30) 


-1 


R = (qT Q) (qT al) 


where the superscripts T and -1 denote the matrix transpose and inverse, 
respectively. 

Equation (30) is unstable in its present form when a small error is added to AL; 
that is, a small error in AL will cause a vastly different solution than when the error 

is zero. This is because the matrix Qi) is almost singular. This condition 
arises because of the high degree of overlap of the weighting functions. Stating this 
another way, the solution is unstable because it depends on information which is not avail- 
able from the measurements. When this solution occurred in the temperature-inversion 
problem, the procedure was to introduce a smoothing matrix and a smoothing parameter 
yw (ref. 32). The value of yw, which depends upon the errors that are present, was 
determined by a trial-and-error procedure. The approach preferred here is to express 
the solution in terms of eigenvectors and eigenvalues. This approach has been used by 
Drayson (ref. 57) and also by Mateer (ref. 58) to deal with the instability of equation (30). 
The advantage of this procedure is that it will show explicitly how many pieces of inde- 
pendent information are available from the measurements. 

The matrix (q t Q) is real, symmetric, and positive definite. Therefore, its eigen- 
values will be real and positive. Suppose that the eigenvalues of (q t Q) designated X^, 
\2, and X 3 are ordered so that Xj = *2 = *3 with corresponding orthonormal eigenvec- 
tors Vj, V 2 , and V3. The column matrix AL) can be expressed in terms of the 
eigenvectors of ( q ? Q) by 


(Q T al ) = ^ b i v i 


(31) 


i=l 


where bj is a constant for each eigenvector. Now, equation (30) becomes 


-1 ° 1 
R = (Q t Q) ^ b^i = ^ bi(Q T Q) Vi 
i=l i=l 


(32) 


It is known by a property of eigenvalues and eigenvectors that 




( 33 ) 
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Therefore, 


3 

B= '>r v i (34) 

i=l xi 


Expanding equation (34) gives 
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(35) 


where ajj gives the components of the eigenvectors v^. That is, 



The errors in the system can be viewed as being contained in the constants bj since, in 
order to determine the values of bj, the matrix AL) is used. The errors arise in Q 
and in AL which, as previously stated, is the difference ^Lmi - Lei), where i is the 
wave-number subscript. Errors in measurement appear through L m j and all other 
errors occur in computing L c i. In reality, therefore, equation (34) can be written 


i=l 1 i=l 1 


(37) 


where b Q i is the value of the constant with no error and e^ represents the error con- 
tribution. The reason for the instability of equation (30) now becomes clear. If Aj is 
very small, then the error term will be very important and could dominate the solution. 
When recalling that the values of A^ are ordered, when equation (37) is expanded, it may 
be assumed that the second and/or third eigenvalues will be so small and the resulting 
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error terms will be so large that the iteration procedure will either be unstable or will 
converge to a physically unrealistic solution. It can be shown (ref. 57) that the eigenvec- 
tors v^ corresponding to these small values of Xj have very little effect upon the cal- 
culated radiance at the top of the atmosphere. Thus, these eigenvectors correspond to 
details of the solution which are easily hidden by a very small amount of noise. Further- 
more, it is clear from equation (37) that inclusion of the small eigenvalues associated 
with these eigenvectors can cause instability or unrealistic results. Therefore, the obvi- 
ous solution is to truncate the terms containing these components. The problem is to 
determine how many terms can be carried in the summation of equation (37) with realis- 
tic errors added to the measured radiances and to various quantities in equation (1). It 
is desirable to carry as many terms as possible since, when a term is dropped, ozone 
information is also lost. There is a trade-off involved between an improved solution due 
to a decrease in noise sensitivity and a degraded result due to loss of information. This 
is discussed further in the next section. The number of terms which can be carried that 
still allow a physically realistic solution is, by definition in this report, the number of 
independent pieces of information available from the measurements. The types and mag- 
nitudes of errors used to determine this information and the responses of the solution to 
these errors are discussed in the subsequent paragraphs. 

ERROR STUDY 

The "measured" radiances for the error study were calculated with use of equa- 
tion (1) by assuming a clear sky and by using vertical ozone and temperature distributions 
measured by the Air Force Cambridge Research Laboratories (refs. 54, and 59 to 61). ^ 
The Planck radiance at the temperature of the ground was used in computing the effective 
lower boundary emission, and the surface emissivity was assumed to be 1.0. Transmit- 
tances used in the computations were determined with the random exponential band model. 
The calculations were limited to the spectral range of 1044 cm - l to 1070 cm - l since this 
is the region where, currently, the absorption-line parameters are known with the great- 
est accuracy. (Refer to fig. 2.) Radiances were calculated at nine wave-number centers 
spaced 3 cm-1 apart. After computing the "measured" radiances, these values were used 
to invert equation (1) to obtain ozone information. This procedure gave complete control 
over all input errors. The "measured" radiances were known perfectly since calculations 
with equation (1) in the inversion procedure were done in exactly the same way as they 

^ These observations were made in situ from a high-altitude balloon at 12 stations 
covering a latitude range extending from the Canal Zone (=9° N) to Greenland (® 77 0 N). 

The ozone measurements were made with chemiluminescent-type ozonesondes developed 
by Regener (refs. 62 and 63). The ozone profile above the maximum balloon altitude is 
determined in this study by extrapolation assuming a constant ozone mixing ratio. 
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were done when the synthetic radiances were originally computed. Convergence criteria 
for these computations were applied to the three perturbation parameters Ap( 03 ) m , 

Ap m , and AH. The iteration was stopped when the following conditions were met simul- 
taneously: Ap(03) m = 1 pmb, Ap m = 1 mb, and AH = 0.01. These limits were selected 
to give sufficient accuracy in a minimum of computing time. If convergence was not 
achieved after nine iterations, the initial guess was altered and the procedure was 
repeated. 

The first series of tests was to evaluate the information content of the radiance 
measurements. Normally distributed random noise was imposed on the "measured" 
radiances, and inversions were performed for various standard deviations ct of the 
noise. The calculations showed that it was not possible to obtain realistic solutions when 
three eigenvectors were included in the expansion of equation (33), even when the radiance 
noise was as low as a = 1 percent. But when the third eigenvector was truncated, the 
inversion was very stable and it gave physically possible solutions for noise as large as 
o = 5 percent. It can be concluded, based on these results, that there are, at most, two 
pieces of independent information available from the measurements. However, it is still 
necessary to show that it is always best to use two terms in the expansion; it may be bet- 
ter to use only one term if the noise level is very high since, as indicated in the analysis 
section, successive truncation of terms continually reduces the sensitivity of the inver- 
sion method to noise. Also, if there are two pieces of independent information available, 
the form of this information must be determined. Both of these tasks were accomplished 
by performing inversions with all input errors equated to zero for various sets of "mea- 
sured" radiances. These calculations were done by using radiances determined for 10 
different ozone profiles covering the latitude range from 9° N to 76.5° N. 

The inversion results showed that in order to obtain physically realizable solutions 
no fewer than two terms should be used in the expansion. Thus, under favorable condi- 
tions, there are two pieces of independent information available from the measurements. 
Furthermore, those two pieces of information are the total ozone u and the altitude h m 
of the maximum in the profile. The data supporting these conclusions are presented in 
table I. Table I shows errors in the total ozone and in the altitude of the maximum of the 
ozone profile when one, two, and three eigenvectors are used in the expansion of equa- 
tion (35). There is very little change in the root-mean-square errors for the 10 profiles 
when two eigenvectors are used rather than three, but there is a drastic change in going 
from two eigenvectors to one. In this latter case, the error in u goes from 0.95 percent 
to 28.8 percent and the error in h m goes from 0.86 km to 3.1 km. It is clear that the 
second eigenvector contains so much information that it cannot be discarded even if the 
noise level is unusually high. The data of table I are explicit examples of the trade-off 
between reduction of noise sensitivity and loss of information when terms of equation (35) 
are truncated. Very little information about u and h m is lost when the third eigen- 
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vector is truncated as shown by the small changes in the root-mean-square errors. Yet, 
the random-noise study showed that when the third eigenvector is included, the inversion 
is so sensitive to noise that the solution becomes unstable even for small errors. There- 
fore, this component should clearly be dropped. On the other hand, a great deal of infor- 
mation is contained in the second eigenvector term and it obviously cannot be neglected. 

The inferred ozone information cannot be expressed explicitly in terms of the three 
parameters of Green's ozone function given by equation (26). The inversion gives one of 
these parameters, p m or h m , but the other piece of information is a function of all 
three parameters through an integral relationship; that is, the total ozone is expressed as 



where p(C>3j is the ozone partial pressure, p is the pressure, and p s is the pres- 
sure at the surface of the earth. 

TABLE I.- ERRORS IN TOTAL OZONE u AND IN ALTITUDE h m OF MAXIMUM 
OZONE PARTIAL PRESSURE WHEN ALL INPUT ERRORS ARE ZERO 


Location 

Date of 
sounding 
used 

Error in u, 
percent, with - 

Error in h m , 
km, with - 

V 1 + v 2 + v 3 

v l + v 2 

V 1 

v i + v 2 + v 3 

v l +v 2 

V 1 

Albrook Fid., 

7/28/65 

-1.3 

-0.68 

58.0 

-0.8 

0.1 

4.8 

Canal Zone 








(9° N) 

11/18/64 

-1.9 

-.78 

44.6 

-.3 

.7 

4.5 

Fla. St. U., 

4/21/65 

-1.4 

-0.72 

12.0 

0.5 

0.6 

1.6 

Florida 








(30.4° N) 

1/11/65 

-1.2 

.61 

37.1 

.4 

.7 

2.8 

Green Bay, 

3/15/65 

0.44 

1.2 

-5.0 

0.75 

0.8 

0.3 

Wisconsin 








(44.5° N) 

2/18/65 

1.2 

-.87 

-7.7 

.6 

.2 

-.5 

Ft. Churchill, 

3/16/65 

0.89 

2.1 

-8.0 

0 

0.2 

-0.2 

Canada 








(58.8° N) 

9/2/64 

-.67 

-.55 

26.2 

.95 

.95 

4.1 

Thule, 

5/19/65 

1.3 

-0.3 

-0.1 

2.8 

2.0 

3.1 

Greenland 








(76.5° N) 

10/7/64 

-.22 

.09 

23.7 

2.3 

.7 

3.8 

Root-mean-square 

1.15 

0.95 

28.8 

1.27 

0.86 

3.1 

errors 
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Since the errors of table I were generated with all input errors equal to zero, these 
errors are the smallest possible values when Green's function is used to model the ozone 
profile. The errors occur because of the mismatch between Green’s function and the 
atmospheric ozone profiles and because of the small amount of information in the mea- 
sured radiances. The actual distribution and the profile determined by inversion using 
two eigenvectors with all errors equated to zero are compared for the July Canal Zone 
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OZONE PARTIAL PRESSURE, plOjKlimb) 


Figure 5.- Comparison of the actual and inferred ozone 
profiles for a July Canal Zone sounding. All input 
errors equal zero. 
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(Albrook Field) sounding and the March Canada (Fort Churchill) sounding in figures 5 
and 6, respectively. These results are typical of solutions obtained for the 10 profiles 
included in this study. The ozone soundings used for figures 5 and 6 were selected 
because they include opposite extremes of latitude, season, and ozone amount. 


TEMPERATURE, T!K) 



Figure 6.- Comparison of the actual and inferred ozone 
profiles for a March Canada sounding. All input 
errors equal zero. 
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The remainder of the error study was done with two eigenvectors in the expansion. 
In the next series of tests, individual errors were imposed one at a time. The types and 
ranges of errors used in these calculations are given in table n. The types of errors in 
table II are those which probably would be encountered in an actual experiment. The 
ranges were selected to facilitate unambiguous identification of the most important error 
effects. The results are presented in figures 7 to 15 for the July Canal Zone profile and 
the March Canada profile. The errors in u and h m from table I have been subtracted 
from the errors in each figure. The only errors plotted for h m are those due to ran- 
dom radiance noise, lower boundary temperature errors, and temperature-profile bias 
errors. The effects on h m resulting from the remaining types of errors were not 
plotted since the deduced values were very small (=0.4 km). 

The error results show a relatively small latitudinal dependence. The errors 
which have the greatest effect on the error in total ozone are random radiance noise, 
lower boundary temperature errors, and absorption-line intensity errors. Those which 
have the greatest effect on h m are random radiance noise, lower boundary temperature 
errors, and temperature-profile bias errors. Fortunately, large temperature bias 
errors are unlikely to occur in practice since the temperature profiles that would be used 
would come from a temperature -inversion method like that of Wark and Fleming (ref. 32) 
which usually gives both positive and negative errors. Also, the other types of errors 
can be controlled. As was pointed out in the analysis section, the effective lower bound- 
ary temperature can be determined to within 0.5 K or less and, undoubtedly, errors in the 
absorption-line intensities will eventually be improved as work advances to eliminate 
defects in current theories regarding the ozone molecule. Once the accuracies of the 
absorption-line intensities are increased, the entire ozone band can be used in the calcu- 
lations thereby causing a reduction in the effective value of the random radiance noise. 


TABLE H.- TYPES AND RANGES OF ERRORS USED IN STUDY OF INDIVIDUAL 
ERROR EFFECTS IN NADIR EXPERIMENT 


Type of error 

Range 

Random radiance noise, normally distributed with a mean of zero, percent . . . 

Standard deviation: 3 o = 1, 3, and 5 
±0.3 x 10- 5 
b 0 ± 20 
±7 
±15 
±2 







a Fifteen inversions were performed for each value of standard deviation. 

b This is the maximum difference between the true profile and the error profile at each of six different alti- 
tudes in the range from the ground to the stratopause. The errors between any two of these altitudes varied linearly 
going from a plus value at one level to a minus value at the next. 
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Figure 9«- Total-ozone error plotted against Figure 10.- Total-ozone error plotted against 

absorption-line intensity error. radiance bias error. 


Total-Ozone Total-Ozone 

Error (%) | Error (%) 




CO KJ- CD 


•ZU T Ft. Churchill, Canada 

(58. 8°N) 

.1.- Total-ozone error plotted against Figure 12.- Total-ozone error plotted against 

temperature profile ± error. temperature -profile bias error. 








Figure 15.- Error in altitude h m of the maximum ozone 
partial pressure plotted against temperature -profile 
bias error. 

To determine the ultimate accuracy which can be expected for u and h m with 

the present inversion scheme, the conditions of table HI were imposed simultaneously. 

The best case values in table HI contain errors which have canceling effects on each 

other, and the worst case values contain errors which are all additive in effect. The 

values represent error limits presently attainable (for the lower boundary temperature 

and temperature -profile errors) or limits which probably can be achieved in the future. 

The standard deviation of the radiance noise was selected in the following way: if the 

standard deviation of the radiance noise is 1 percent, which is currently possible, and 

if 30 wave numbers can eventually be used (once the accuracies of the absorption-line 

intensities are improved), then the effective standard deviation can be reduced to 

ff = -7= = - ‘0 - = 0.18 percent. This value can be simulated in the current calculations 
e i/n \/30 

with o = 0.54 percent, since with nine wave numbers <7 e ff = ^^4 = 0.I8 percent. The 



TABLE m.- ERRORS USED IN SIMULTANEOUS ERROR STUDY OF NADIR EXPERIMENT 


Type of error 

Best case values 

Worst case values 

Radiance noise, percent 

a = 0.54 

a = 0.54 

Radiance bias, watts/m -sr 

-0.15 x 10“ 5 

0.15 x 10-5 

Temperature profile, alternating sign, kelvins . . 

±3 

±3 

Ozone absorption-line intensity, percent 

4 

4 

Lower boundary temperature, kelvins 

0.4 

-0.4 


TABLE IV.- ERRORS IN TOTAL OZONE u AND IN ALTITUDE h m OF MAXIMUM 
OF OZONE PROFILE WHEN ERRORS OF TABLE HI 
ARE IMPOSED SIMULTANEOUSLY 


Location 
and date of 
sounding used 

Best case errors 

.t 

Worst case errors 

Error in 
u, percent 

Error in 
h m , km 

Error in 
u, percent 

Error in 
h m , km 

Mean 

Std. 

dev. 

Mean 

Std. 

dev. 

Mean 

Std. 

dev. 

Mean 

Std. 

dev. 

Albrook Fid., 
Canal Zone 
(9° N) 
7/28/65 


■ 

0.17 

0.3 

-16.0 

3.5 

-1.0 

0.4 

Ft. Churchill, 
Canada 
(58.8° N) 
3/16/65 

■ 

3.4 

0.15 

0.16 

-14.3 

3.4 

0.28 

0.2 


remaining quantities of table m are estimates of final accuracies. The resulting errors 
are given in table IV. The mean error in total ozone ranges from 3.4 percent in the best 
case to -16 percent in the worst case with a standard deviation of about 4 percent. Thus, 
it is not unreasonable to expect that total ozone can eventually be determined to within 
10 percent or less. Also, the error in h m varies from 0.15 km in the best case to 
-1 km in the worst case with a standard deviation of approximately 0.3 km. It would 
appear that this parameter can ultimately be determined to within about 0.5 km. How- 
ever, some caution should be exercised here since table I shows several errors in h m 
which are greater than 0.5 km even when all input errors are zero. The "all error" 
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study indicates such small discrepancies because the two ozone profiles used in the study 
were matched reasonably well by Green's ozone function (eq. (26)), thereby giving a small 
"zero error" for h m in both cases. This will not always be true as table I shows. How- 
ever, it is reasonable to expect, based on tables I and IV, that h m can be determined to 
within about 1.5 km, at least for the large majority of profiles encountered in the 
atmosphere. 

The calculations of table IV were repeated for varying degrees of cloudiness in the 
troposphere. To simulate the effect of clouds, the effective lower boundary emission was 
computed by weighting emissions from the ground and cloud tops according to the per- 
centage cloud cover. 4 The clouds were considered to be black, and the Planck radiance 
B e ff of the effective lower boundary was computed according to the equation 

Beff = Bgnd(l - Pql) + ®cL PqL (39) 

where Bg n( j is the Planck radiance at the ground temperature, B C L is the Planck 
radiance at the cloud-top temperature, and P c l is the percentage cloud cover. Alto- 
stratus cloud covers of 100 percent and 8 percent at 5 km were used for the Canal Zone, 
and values of 100 percent and 12 percent at 3 km were used for the Canada calculations. 
These conditions, exclusive of the complete cloud covers, are averages for the respective 
locations according to data published by London (ref. 64). The calculations revealed only 
a small effect due to the presence of clouds. Negligible changes occurred in the errors 
of table IV when partially cloudy conditions were used, and the errors increased only 
slightly when a complete cloud cover was used. The means and standard deviations of the 
errors in u increased by about 2 to 3 percent, and the errors in h m increased about 
0.2 km. Therefore, it appears that the presence of clouds will not seriously degrade the 
inversion results as long as there is some contrast between the ozone spectrum and the 
effective lower boundary emission spectrum. (See appendix B.) 

Since the band model used in the error study did not allow for variation in 
absorption-line intensity with temperature, a test was conducted to evaluate this effect. 

The line-by-line integration technique was used to compute the atmospheric transmit- 
tances for a mid-latitude ozone profile and a corresponding temperature profile. This was 
done for the reference temperature profile and for a profile 10 K higher in temperature at 
every point. These two sets of transmittances were then used along with the reference 
temperature profile and equation (1) to compute the radiances at the top of the atmosphere. 
The results are plotted in figure 16. The computations show a slight temperature depen- 
dence, mostly in the wings of the ozone band, but the effect is so small that it would not 
change the results of the error study significantly if it were included in the calculations. 

4 The percentage cloud cover is not required in applying the inversion technique to 
flight data since cloud effects are treated indirectly by the method discussed in the anal- 
ysis section. 
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Figure 16.- Effect of temperature on the ozone spectrum. 

Dashed curve computed by using reference temperature 
profile T r (p) . Solid curve computed by using refer- 
ence temperature profile plus 10 kelvins. 

The inferred values of u and h m would differ by about 8 percent and 0.4 km, respec- 
tively, for the two spectra of figure 16. This effect is large enough that it should not be 
neglected in precise calculations such as when applying the inversion method to actual 
radiances measured from a balloon or satellite. 

APPLICATIONS 
Palestine, Texas, Data 

The analysis methods were applied first to radiances measured from a balloon at 
the 31-km level above Palestine, Texas (31.8° N, 96° W). The measurements were made 
on May 8, 1966, in a flight conducted by the High Altitude Engineering Laboratory of the 
University of Michigan. The instrumentation and other details concerning the flight have 
been discussed in a report by Chaney et al. (ref. 29). A breadboard version of an IRIS 
having a resolution of 5 cm~l was used to make the radiance measurements, and a typical 
spectrum from the experiment was presented in figure 1. 




The data were of very good quality. The standard deviation of the radiance noise 
was only 0.5 x 10"^ watt/m-sr which corresponds to a level of about 1.3 percent in the 
center of the ozone band and less than 1 percent in the wing region. This low noise made 
\ the data especially well suited to test the inversion method. Unfortunately, neither the 
total ozone nor the vertical ozone distribution was measured at the time and location of 
the balloon. Therefore, it was not possible to make a direct comparison between mea- 
sured quantities and values determined from the inversion. However, the total ozone 
could be estimated for Palestine based on a Dobson measurement made at Albuquerque, 
New Mexico (35° N, 106.6° W). According to H. S. Muench of the Air Force Cambridge 
Research Laboratories, the ozone amount at Albuquerque on the day of the balloon flight 
was 0.317 atm-cm, and since Palestine, Texas, is located at a slightly lower latitude, 
the total ozone there should have been slightly less; a reasonable guess would be a 
value of 0.298 atm-cm. This estimate is based on a total-ozone latitudinal gradient 
of 0.006 atm-cm per degree of latitude as determined from data presented by Hering 
and Borden (ref. 61). There was no vertical sounding at Albuquerque on the day of the 
balloon flight which could be used to estimate the altitude of the primary maximum in 
the ozone profile. Past measurements at the same time of year at Albuquerque indi- 
cate that the value of h m could lie anywhere in the range from about 23 km to 27 km. 

Transmittances used in the inversions were calculated by the line-by-line method 
rather than by the less accurate random-exponential-band-model approach. Inversions 
were attempted by using the radiances of table V and the atmospheric temperature profile 
measured in situ from the balloon. A stable solution could be obtained only after water- 
vapor line absorption was neglected in determining the lower boundary emission. This 
was equivalent to using background radiances which were too high; that is, it was equiv- 
alent to having a positive error in the lower boundary temperature. The most likely rea- 
son that the higher background emission was needed was that it compensated for a positive 
absorption-line intensity error. Figures 8 and 9 show that positive errors in the lower 
boundary temperature and the absorption-line intensity have canceling effects on each 
other. A positive absorption-line intensity error should be expected based on figure 2 
which shows that there is too much calculated absorption over the spectral range of 
table V for a homogeneous path. Other calculations with the Nimbus data also show that 
this is true for an atmospheric path. 

The theory of canceling error effects is plausible when estimates of the errors 
involved are considered. For example, the background radiance error was about 
lO" 5 watt/m-sr in the wings of the band. This corresponds to a lower boundary temper- 
ature error of about 1 K. According to figure 8, a 1 K error, causes a total-ozone error 
of about 25 percent. The error in absorption-line intensity can be estimated from fig- 
ure 1. A representative error for the interval of table V is about 35 percent. If the 
curve of figure 9 is extrapolated, it shows that a 35-percent error causes a total-ozone 
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error of around -30 percent. When the two sources of error are combined, the net result 
is a total -ozone error of only about -5 percent. 


TABLE V.- PALESTINE, TEXAS, DATA 


Wave number, cm-l 

Radiance, L(£,0), watts/m- sr 

1052 

37.2 x 10" 5 

1056 

38.6 

1060 

45.2 

1062 

50.3 

1066 

61.5 

1068 

66.8 

1070 

70.8 


There was almost complete cancellation of error effects for the spectral range of 
1060 cm-l to 1070 cm~l where the absorption-line intensity error was large and positive. 
However, this was not true for the remaining region of table V (1052 cm-l to 1060 cm-l) 
where the absorption-line intensity errors were different. Figure 2 shows that the 
absorption errors are smaller in the range of 1052 cm-l to 1060 cm-l and the error 
changes sign; at 1052 cm-l the error is almost zero and it is negative. Thus, when the 
wave numbers centered at 1052 cm-l and 1056 cm-l were included in the calculations, 
the balance of error effects was upset. The effect of using these wave numbers was to 
cause the computed spectrum to shift upward in the wing region of the band. The shift 
was due to a change in slope of the spectrum in the wing region brought about by the least- 
squares fitting procedure. The slope became progressively steeper as first 1056 cm-l 
and then 1052 cm-l was included. The steeper slope gave higher radiances at every point 
which corresponded to less ozone absorption and, hence, smaller ozone amounts. For 
example, when all wave numbers of table V were used, the inverted ozone amount was 
0.129 atm-cm; when the 1052 cm-l wave number was dropped, the value of inferred total 
ozone increased to 0.161 atm-cm. Finally, when both wave numbers were dropped to 
give good cancellation of error effects, the value of inferred total ozone increased to 
0.294 atm-cm. Since the balloon was not at the top of the atmosphere, this value should 
be increased by about 5 percent to account for the ozone above the balloon. When this is 
done, the inferred total ozone becomes 0.309 atm-cm. This is in excellent agreement 
with the value estimated for Palestine, Texas. The Green ozone function (eq. (24)) corre- 
sponding to this inversion is shown in figure 17. The value of h m is 26.6 km which is 
within the expected range of values for the latitude and time of year of the measurement. 
These results, coupled with the results of the error study, prove the validity of the inver- 
sion method developed in the analysis section. 
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OZONE PARTIAL PRESSURE, plOjHumb) 


Figure 17. - Ozone profile inferred by using radiances measured from 
a balloon above Palestine, Texas, on May 8, 1966. 

The results also suggest a need for further work to increase the accuracy of the 
absorption-line intensities. Because of the low noise level in the Palestine radiances, it 
was possible ultimately to perform the inversion by using only five wave-number centers 
in a narrow spectral interval (1060 cm-1 to 1070 cm-1). With this narrow interval, a 
cancellation of error effects due to absorption-line intensity errors was obtained by add- 
ing positive background radiance errors. However, such a narrow interval will not afford 
enough suppression of random radiance noise to permit accurate results for noise levels 
normally encountered (this is vividly illustrated in later calculations with the Nimbus HI 
data); and when a wider interval is used to reduce the effective noise level, error effects 
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no longer cancel because of the change in algebraic sign of the absorption error in one- 
half of the band. This illustrates the need for more accurate values of the absorption 
parameters across the band. 


Nimbus Data 

The solution procedures were also applied to radiances measured by an IRIS 
mounted on the Nimbus HI satellite launched in April 1969. These data were kindly pro- 
vided to the author by B. J. Conrath of the Goddard Space Flight Center. The location of 
the subsatellite point at the time of the measurements was 20.8° N, 70.9° W. This is 
very close to the Grand Turk Air Force Base (21.5° N, 71.1° W) where an ozonesonde 
observation was made less than 24 hours before the satellite pass. Since daily ozone 
variations are quite small in the temperate zone, the ozone profile probably changed only 
slightly between the time of the ozonesonde launch and the time that the radiances were 
measured. Therefore, the ozonesonde data should be a good gage for evaluating the inver- 
sion results. 

Unfortunately, the Nimbus HI measurements were of poor quality in the vicinity of 
1042 cm - l. The noise-equivalent radiance was about 2 x 10"5 watt/m-sr (ref. 28). 

This corresponds to an error of about 5 percent in the center of the ozone band and of 
approximately 3 percent in the wing region. Considering the mean error to be 4 percent 
gives a 3a radiance error of 12 percent; and, according to figures 7 and 13, this causes a 
3a total-ozone error of about 96 percent and a 3a h m error of about 7.7 km. There- 
fore, it is anticipated that inversion results using the Nimbus data will be poor. 

Inversions were performed by using the data of table VI and the temperature profile 
measured in situ at the time of the ozonesonde observations. The first inversion was 
attempted by using background radiances that included water-vapor line absorption and, 
as with the Palestine data, the solution was unstable. Therefore, water -vapor line 
absorption was neglected for the same reason as was discussed in the previous subsection. 

TABLE VI.- NIMBUS DATA 


Wave number, cm-l 

Radiance, L(£,0), watts/m -sr 

1053.3 

40.3 x 10- 5 

1055.4 

42.1 

1059.5 

45.4 

1061.6 

49.0 

1065.8 

56.4 

1067.9 

59.5 

1070.0 

60.2 
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When wave numbers less than a value of about 1060 cm-1 were used, the solution was al/30 
unstable. This was probably due to the combined influence of the high-radiance noise and 
the imbalance of error effects brought about by absorption-error differences in the spec- 
tral region below 1060 cm-1 (as discussed in the previous subsection). Finally, an inver- 
sion was performed by using the same spectral range as was ultimately used with the 
Palestine data and the solution was barely stable. It did not converge until the seventh 
iteration. Usually, convergence is achieved in only three or four iterations. The inver- 
sion gave a total-ozone value of 0.062 atm-cm and a value for h m of 18.7 km. These 
correspond to errors of -76 percent and -6.3 km, respectively, which, as expected, are 
large errors. The error magnitudes are consistent with the values predicted from fig- 
ures 7 and 13. These results could be improved by using the entire ozone band to reduce 
the effective noise level, but as previously discussed, existing absorption errors prohibit 
this. 

In order to evaluate typical absorption errors for an atmospheric path, the radiance 
at the top of the atmosphere was computed by using ozonesonde and temperature data 
taken at Grand Turk. Water-vapor line absorption was considered in computing the back- 
ground radiances, and the line-by-line technique was used to compute the transmittances. 
Computed values were compared to the radiances measured at the satellite and the results 
are shown in table VII. The quantity AL = L m i - L C i is the difference between the mea- 
sured and calculated radiances. Table VII shows that there is too much calculated 
absorption in the region of 1060 cm-1 to 1070 cm* 1 for the atmospheric path just as it is 
for the homogeneous path spectrum shown in figure 2. This supports the argument of the 
previous subsection that justified the need for positive background radiance errors 
because of the presence of positive absorption-line intensity errors. 

TABLE VII.- DIFFERENCES BETWEEN RADIANCES MEASURED ON-BOARD NIMBUS HI 
AND RADIANCES COMPUTED BY USING CORRESPONDING OZONE PROFILE 


Wave number, cm -1 

Radiance difference, AL, watts/m-sr 

1053.3 

-1.6 x 10" 5 

1055.44 

.33 

1059.5 

1.5 

1061.6 

4.1 

1065.8 

4.7 

1070 

1.4 
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CONCLUDING REMARKS 


A workable technique has been developed for inverting nadir measurements to obtain 
ozone information. The success of the method depends critically on the accuracy of the 
quantity which represents the energy emitted from the lower boundary under a satellite. 

A method was formulated for computing this quantity which circumvents the troublesome 
problems caused by the water-vapor continuum absorption and by atmospheric clouds. 

By incorporating water-vapor line absorption in the calculations, the effective tempera- 
ture of the lower boundary could be found to within 0.5 K or less. The inversion equations 
were expanded in terms of the eigenvectors and eigenvalues of a least-squares -solution 
matrix, and an analysis was performed by using synthetic clear-sky radiances computed 
for 10 different ozone profiles. The results showed that, under favorable conditions, 
infrared radiance measurements in the ozone band contained two pieces of independent 
information: the total ozone u and the altitude h m of the primary maximum in the 
ozone profile. It is not possible to obtain the value of the maximum ozone partial pres- 
sure or the width of the ozone profile from these data. A comprehensive study, using a 
variety of error sources, revealed that errors in u are most sensitive to random radi- 
ance errors, lower boundary temperature errors, and ozone absorption-line intensity 
errors. Errors in h m are affected most by the former two errors and by temperature- 
profile bias errors. 

In the next series of tests, all error sources were considered simultaneously. The 
results showed that it should ultimately be possible to determine u to within 10 percent 
or less and to determine h m to within 1.5 km when the root-mean-square radiance noise 
level is 1 percent or less. These calculations were also made for varying degrees of 
cloudiness in the troposphere. The data showed that the presence of clouds does not 
seriously affect results as long as there is some contrast between the ozone spectrum 
and the effective lower boundary emission spectrum. Sometimes, in high-latitude areas 
of the globe, there may be little or no contrast between these spectra due to the effects of 
temperature profile and cloud conditions. Under these circumstances, it may be difficult 
or impossible to infer ozone information from the measurements. 

Finally, the inversion technique was applied to balloon data measured over 
Palestine, Texas, and to the Nimbus HE satellite data measured over Grand Turk Island 
in the Bahamas. The two sets of measurements were in sharp contrast with each other; 
the balloon data contained very little noise, whereas the satellite data were very noisy by 
comparison. The inversion results reflected this contrast. The Palestine inversion was 
satisfactory considering the uncertainties which existed. The errors in u and h m 
were estimated to be well within the limits of 10 percent and 1.5 km, respectively. How-' 
ever, the Nimbus results were poor. The total ozone u was off by 76 percent and h m 
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was in error by 6.3 km. If absorption errors throughout the band had been smaller, the 
Nimbus results could have been substantially improved by using the entire ozone band to 
reduce the effective radiance noise level. Unless absorption-line intensities are 
improved by empirical or other methods, calculations will have to be limited to a narrow 
spectral interval and, consequently, to low-noise data (less than 1 percent). 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 13, 1973. 
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APPENDIX A 


RELATIONSHIPS FOR THE RANDOM EXPONENTIAL BAND MODEL 
AND THE CURTIS-GODSON APPROXIMATION 


The transmittance through a homogeneous layer is given by Goody (ref. 42) as 


= 


exp 


r 





1 + 2 






(Al) 


and he shows that, in order to force a fit in the weak-line (u^. « 1 
(u k » regions, the following expressions should be used: 


and strong-line 


d 



-> 2 


H 

2 ^Lo^) 



(A2) 


(A3) 


where N' is the total number of lines in the interval Av and o^ 0 is the Lorentz half- 
width of the ith line at standard temperature and pressure. 

According to the Curtis-Godson approximation the effective pressure of the atmo- 
sphere between the satellite and an arbitrary pressure p Q is 


Pk = 



[t(p) 

If* 

JN 

T(pj] 

f* 


(A4) 
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APPENDIX A - Concluded 


where q is the absorber mass mixing ratio and g is the acceleration of gravity. The 
effective absorber amount is 


u k = 



s[T(pj)|dp 

W) 


(A5) 


where T is a mean temperature. However, if s|t(p) 1 = S^T^, that is, if the line inten- 
sity does not vary with altitude, then equations (A4) and (A5) become 


Pk = 


u k = 



(A6) 


(A7) 


The effective length of the homogeneous path is computed by 



(A8) 


These last three equations, along with equations (Al), (A2), and (A3), are used throughout 
this report to compute the transmittance when the random -exponential band model is used. 
The line intensities required to compute equations (A2) and (A3) were taken from Clough 
and Kneizys (ref. 35). The half-width was considered to be the same for all lines and 
equal to 0.08 cm~l at standard temperature and pressure. This value was suggested by 
Clough and Kneizys (ref. 35) and used by Goldman and Kyle (ref. 41). 
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APPENDIX B 


EFFECT OF LOWER BOUNDARY EMISSION ON THE OZONE SPECTRUM 


The energy emitted by the underlying boundary can have a significant effect on the 
ozone spectrum measured at the top of the atmosphere. Assume, for the present, that 
the emissivity of the lower boundary is one and that the ozone layer is isothermal at a 
temperature T(p r ^ which is different from the temperature T^p^) of the lower bound- 
ary where p r is the pressure at the bottom of the ozone layer and Pg is the pressure 
at the lower boundary. These conditions are represented in figure 18. In this case equa- 
tion (1) can be written as 


1^0,0) = B i>o,T(p g ) T^ 0 , Pg ) - B v 0 ,T(p r ) ^ °’ P ^ d^v 0 ,p) 

r^o.Pg) 

^oPt) 


v 0 ,T(p) dr(i> 0 ,p) 


(Bl) 


where prj. is the pressure at the top of the effective lower boundary. If it is assumed 
that t^ 0 ,Pt) ~ T (i'o>Pg)> then equation (Bl) reduces to the following after integration: 



(B2) 


The spectrum defined by equation (B2) and different values of v 0 will show the structure 
of the ozone band because of the factor rf^pi However, if the lower boundary is at 
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APPENDIX B - Concluded 


the same temperature as the ozone layer, that is, if T^Pg') = T^p r ^, then equation (B2) 
reduces to 


L(>0) = B 5> 0 , T (pr) 


(B3) 


which shows no structure. Such a spectrum contains no ozone information. Ways in 
which equation (B3) can occur are illustrated in figures 19 and 20. In figure 19 the under- 
lying boundary is the earth's surface, and in figure 20 it is the top of a cloud layer for a 
completely overcast condition. If the troposphere is partly cloudy the effective lower 
boundary temperature, for a temperature profile like that of figure 20, will be at some 
intermediate value between the ground and tropopause temperatures. There would be 
structure in the ozone band for this partly cloudy condition, but the structure would not 
be as pronounced as for the clear-sky case. In fact, equation (B2) shows that as the 
value of T^Pg!) approaches the value of T^p r ^ the structure in the ozone band decreases 
and finally disappears when the two temperatures are equal. 

A broad, nearly isothermal region in the middle and lower stratosphere is typical 
of the temperature profile in the high-latitude areas of the globe. Therefore, the entire 
range from full to no structure in the ozone band may be encountered at high latitudes 
depending upon the temperature profile and cloud conditions. Consequently, it sometimes 
may be difficult or impossible to determine ozone information from satellite measure- 
ments in the high -latitude regions. 




TEMPERATURE 


Figure 19.- Isothermal ozone layer with 
T (p r ) = T(p g ). 


Figure 20.- Isothermal ozone layer with 
a total cloud cover in the tropo- 
sphere. T^p r j = T (Pg)- 
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